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We present a coupled system of integral equations for the KK — > NN and KK — > NN S-waves 
derived from Roy-Steiner equations for pion-nucleon scattering. We discuss the solution of the 
corresponding two-channel Muskhelishvili-Omnes problem and apply the results to a dispersive 
analysis of the scalar form factor of the nucleon fully including KK intermediate states. In par- 
ticular, we determine the corrections A a and A^, which are needed for the extraction of the pion- 
nucleon a term from nN scattering, and show that the difference Ad — Ao- = (— 1.8±0.2)MeV 
is insensitive to the input nN parameters. 
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1. Introduction 

The pion-nucleon a term c k n measures the contribution of the light quarks to the nucleon 
mass m, and is directly related to the form factor of the scalar current 

o(t) = ^-(N(p')\m(uu + d-d)\N(p)), t = {p'-p)\ m='^±^, (1.1) 

at vanishing momentum transfer a(0) = g k n. The standard procedure for its extraction from pion- 
nucleon (nN) scattering relies on the low-energy theorem [[T], ^] 

F*D+(s = m 2 ,t = 2Ml) = a(2M 2 n ) + A R , (1.2) 

which relates the Born-term-subtracted isoscalar nN scattering amplitude at the Cheng-Dashen 
point D + (s = m 2 ,t = 2M 2 ) to the scalar form factor evaluated at 2M 2 . The remainder Ar is free 
of chiral logarithms at full one-loop order in chiral perturbation theory (ChPT) [|],^J], and has been 
estimated as ^ 

|A R |<2MeV. (1.3) 



Rewriting (1.2) in terms of 



A D = F 2 {D + (s = m 2 ,t = 2Ml)-d+-2M 2 K d+y A a = o(2M 2 n ) - a nN , (1.4) 

the extraction of the a term reduces to the determination of the subthreshold parameters d^ and 
as well as the combination Ad — A a — Ar. The first two corrections can be calculated using a 
dispersive approach [§] 

A D -A a = (-3.3±0.2)MeV, (1.5) 

where the error only covers the uncertainties in the nn phase shifts available at that time. Here, we 
update the determination of A^ and A a using modern nn phases, fully including KK intermediate 
states, and carefully studying the dependence of the results on nN subthreshold parameters as well 
as the 7T./V coupling constant. 



2. Scalar pion and kaon form factors 

We first consider the case of the scalar pion and kaon form factors F%(t) and F^(t), which 
serve both to illustrate the method and as input for the scalar form factor of the nucleon. Unitarity 
in the nn/KK system intertwines both form factors according to ^ 

lmF s (t) = (T(t))*L(t)F s (tl F s (t) = , CD 

with the phase-space factor 

Z{t) = diag(c?e(t-t M ),(j! c O(t-t K )y a}=Jl-j, h=AM 2 ie{n,K}, (2.2) 
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Figure 1: Modulus and phase of the scalar pion form factor. The solid, dashed, and dot-dashed lines refer to 
f£(0) = M\j1, 0.4 M\, and 0.6 M\. The phase of F%(t) is compared to 5°, as indicated by the dotted line. 



and the T-matrix 



Tit) 



(t)\e^ 



(2.3) 



2ia, K 



expressed in terms of the KK and kk 



KK phase shifts 8$ and \f/$ as well as the inelasticity pa- 



rameter T]° = J\ -4a, 7r a/ r |g(f)| 2 0(? -t n ). The two-channel Muskhelishvili-Omnes (MO) prob- 
lem [|7|, |8]| defined by the unitarity relation ( |2.1| ) permits two linearly independent solutions Qi, 
^2 [Q], which may be combined in the Omnes matrix £l(t). In general, there is no analytical 
solution for £l(t), we follow here the discretization method of [Q] for its numerical calculation. 

Since the form factors are devoid of a left-hand cut, they are related directly to the solutions of 
the MO problem with coefficients determined by F%(0) and F|(0) [0]. Using ChPT at @{p A ) and 



the low-energy constants from Q1Q[ ] we find 



F*(0) = (0.984 ±0.006)M^ f|(0) = (0.4. . .0.6)m£, 



(2.4) 



which, together with 5q an d Ho fr° m an extended Roy-equation analysis of KK scattering [11], 
from partial- wave analyses [ 12 ] , and \g(t)\ from a Roy-Steiner (RS) analysis of kK scattering []l3|] , 
yield the results for F%(t) depicted in Fig. [jj The strong dependence of F%(t) near t% on F^(0) at- 
tests to the inherent two-channel nature of the problem and implies that an effective single-channel 
description in terms of the phase of F%(t) only works for sufficiently large F^(0). 



3. From Roy-Steiner equations to the scalar form factor 



Unitarity couples the KK —?- NN and KK —?■ NN 5-waves /V (?) and lr + (i) analogously to ( |Ll| ) 



hnm = {T(t)Ymm, m 



(3.1) 



but due to the presence of the left-hand cut the solution of the corresponding MO problem involves 
inhomogeneous contributions, which may be derived from RS equations, cf. [[T3|-[l5|]. Generically, 
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Figure 2: Results for the real and imaginary part of f%(t). The solid, dashed, and dot-dashed lines refer to 
the input RSI, RS2, and RS3 as described in the main text. The black crosses indicate the results of [|l7|]. 



the integral equation takes the form 

f(t)=A(t) + (a + bt)(t-Am 2 ) + 



2 t 2 (t-Am 2 ) 



71 



dt' 



Im/(0 



t' 2 (t' -4m 2 )(t' -t)' 



(3.2) 



where A(f) includes Born terms, ^-channel integrals, and higher ^-channel partial waves, while a 
and b subsume subthreshold parameters that emerge as subtraction constants. The main difficulty 
in the evaluation of the formal solution 



fit) = A(t) + (t- Am 2 )Q.(t) (1 -tQ.(0))a + t{t- Am 2 )Q.(t)b 



(3.3) 



t 2 (t 



71 



Am 2 ) „ f . ImQ-Ht'W) 

-a(t) / dt'- 



t' 2 (t'-Am 2 )(t'-t 



t 2 (t-Am 2 ) „ 
+ — J -Slt) 



71 



dt 



, a-\t')\mf{t') 

t' 2 (t' -Am 2 )(t' -tY 



concerns the construction of the Omnes matrix for a finite matching point t m jr4|]. 

In the numerical analysis we put Im/(?) = above t m , which we choose as t m = Am 2 (thus 
exploiting a kinematical zero of f(t)), take the 71^ and KN ^-channel partial waves from [jl|], 
and use the KH80 7cN coupling constant and subthreshold parameters as reference point [|n|]. 
In order to assess the uncertainties for higher energies we consider the following variants of the 
input: first, we keep the phase shifts 8$ and i//q constant above ^/Iq = 1.3 GeV ("RSI"), where An 
intermediate states become important and thus the two-channel approximation will break down, 
and second, guide both phase shifts smoothly to their asymptotic value of 2% as for the meson 
form factors ("RS2"). Finally, we amend RSI in such a way that ^{t), the KN component of 
the inhomogeneity, is put to zero in order to assess the uncertainty in the K/V input ("RS3"). The 
corresponding results for /+(?) depicted in Fig. || show that the largest uncertainty is induced by 
the high-energy phase shifts. 



4. Results 

The scalar form factor of the nucleon fulfills the unitarity relation 
Ima(f) 



(4.1) 
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Figure 3: g k n an d A a as a function of the integral cutoff A. 



so that, based on the results of the previous sections, the un- and once-subtracted dispersion rela- 
tions 



alt) = - / dt' , ' =o KN + - / At' - , 

W 71 J t'-t 71 J t'(t'-t) 



(4.2) 



tit r. 



evaluated at t = and t = 2M 2 in principle determine o k n and A a provided the two-channel ap- 
proximation for the spectral function is sufficiently accurate in the energy range dominating the 
dispersive integral. Fig. ||| shows that, while the dispersion relation converges too slowly for the 
a term itself, the result for A a becomes stable for A > 1 GeV. Adding the uncertainties from the 
spectral function and the variation of the integral cutoff between A = 1.3 GeV and A = 2m, we find 

A a = (13.9 ± 0.3) MeV 

7+ il/f I 1 AA \ I 7. I A+ fl/f3 1 1/11 i 7 . I U+ fl/f3 



+ZA 14.28 1 +Z 2 ^+M W + 1.46J +Z 3 [d+Mj c -l.U) + Z 4 (£+ M^ + 3.54 
Zi = 0.36MeV, Z 2 = 0.57MeV, Z 3 = 12.0MeV, Z 4 = -0.81MeV, (4.3) 

where we have made the dependence on the 7iN parameters explicit (note that more modern deter- 



minations point to lower values of the 7iN coupling constant around g 2 /4n ~ 13.7 [18-20]). 
Similarly, the correction Ap follows from the ^-channel expansion 



f lm f (t') 

D + (s = m 2 ,t) = 4 +4l' - \6t 2 J At ' jri^, — 4^2)^/ — 7J + i J - 2 ) + {^-channel integrals} 

(4.4) 



evaluated at t = 2M 2 ., which gives 

A D = (12.1 ±0.3) MeV 

+Zj - 14.28^ +Z 2 (d+M n + 1.46) +Z 3 (4*4 - I.I4) +Z 4 (b+Ml + 3.54 
Zi = 0.42MeV, Z 2 = 0.67MeV, Z 3 = 12.0MeV, Z 4 = -0.77MeV. (4.5) 
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Comparison with ( |Q| ) shows that the dependence on the 7iN parameters cancels nearly completely 
in the difference 

A fl -A ff = (-1.8±0.2)MeV. (4.6) 

This cancellation can be explained by the observation that the spectral function in both dispersion 
relations involves /? in a very similar manner, so that both integrals are equally affected by the 
dependence on the 71N parameters inherited from f^_. In the same way, part of the uncertainties 
discussed in Sect. |3| drop out, so that the final error estimate for Aq — A a even decreases compared 
to the uncertainty in both corrections individually. 
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